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Transmitting  antennas  located  near  the  earth  have  a  problem 
with  electromagnetic  scattering  caused  by  the  ground.  Researchers 
suggest  that  the  inverse  of  this  problem  may  be  used  to  determine 
the  complex  permittivity  of  the  ground.  That  application  is  the 
goal  of  this  thesis. 

The  solution  of  the  half-space  problem  was  originally  formula¬ 
ted  in  1909  by  Arnold  Sommerfeld.  Unfortunately,  his  formula¬ 
tion  was  written  in  terms  of  a  set  of  integrals  that  tend  to  have  a 
slow  rate  of  convergence.  Since  then,  researchers  have  sought 
faster  methods  of  computing  fields  in  the  presence  of  an  arbitrary 
infinite  half  space  by  making  approximations  to  Sommerfeld's 
integrals  to  increase  their  convergence  rate.  Others  have  replaced 
the  ground  with  an  equivalent  image  current. 

Much  of  the  work  on  this  subject  has  been  directed  toward  the 
inverse  scattering  problem,  i.e.,  inferring  the  dielectric  properties 
of  the  ground  by  measuring  scattered  fields,  using  time-domain 
techniques  or  continuous-wave  measurements. 

In  choosing  one  method  over  another,  the  parameters  related 
to  the  properties  of  the  ground  should  be  large  enough  or  small 
enough  to  be  measured  with  reasonable  accuracy.  The  experi¬ 
mental  setup  should  be  straightforward  and  simple,  and  the 
analysis  should  not  be  overly  complex. 

Relating  the  input  impedance  of  a  dipole  near  the  ground  to 
its  permittivity  was  deemed  a  practical  experimental  technique. 
The  input  impedance  of  a  half-wave  dipole  is  around  75  ohms, 
well  within  the  measurement  range  of  modern  RF  impedance 
meters,  and  the  experiment  requires  only  a  one-port  measurement 
of  a  simple  antenna.  The  analysis  involved  was  manageable.  The 
method  chosen  for  theoretical  analysis  was  an  image  approach 
that  replaces  the  ground  with  an  equivalent  image  current  located 
on  a  line  in  complex  space. 

This  study  addresses  the  theory  between  the  electric  properties 
of  the  ground  to  the  input  impedance  of  a  horizontal  electric  di¬ 
pole,  including  numerical  computations,  patterns  of  the  input 
impedance,  practical  considerations,  and  recommendations  for 
experimental  methods  of  determining  the  ground  permittivity 
using  this  theory. 
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1.  INTRODUCTION 


The  problem  of  electromagnetic  scattering  caused  by  the  ground  has  been 
investigated  by  a  great  number  of  researchers  over  the  past  few  decades.  It  is  an  important 
problem  because  transmitting  antennas  are  often  located  near  the  earth.  Many  have 
suggested  that  the  inverse  of  this  problem  may  be  used  as  a  means  of  determining  the 
complex  permittivity  of  the  ground  [1,2, 3, 4 ,5].  The  latter  application  is  the  goal  of  this 
paper. 

The  solution  of  the  half-space  problem  was  originally  formulated  in  1909  by  Arnold 
Sommerfeld  [6].  His  formulation  is  written  in  terms  of  a  set  of  integrals  now  commonly 
known  as  Sommerfeld  integrals.  Unfortunately,  these  integrals  tend  to  have  a  slow  rate  of 
convergence.  Since  that  time,  many  researchers  have  sought  faster  methods  of  computing 
fields  in  the  presence  of  an  arbitrary  infinite  half  space.  Some  have  made  approximations 
to  Sommerfeld's  integrals  to  increase  their  rate  of  convergence  [7,8,9].  Others  have 
replaced  the  ground  with  an  equivalent  image  current  [10,1 1,12,13]. 

Many  of  the  papers  on  this  subject  have  been  specifically  directed  toward  the 
inverse  scattering  problem,  i.e.,  inferring  the  dielectric  properties  of  the  ground  by 
measuring  scattered  fields.  Some  suggested  time-domain  techniques  [3,5].  Others 
discussed  continuous- wave  measurements  [1,2,4]. 

In  choosing  one  method  over  another,  a  number  of  considerations  need  to  be  taken 
into  account.  One,  the  parameters  related  to  the  properties  of  the  ground  should  be  easy  to 
measure,  i.e.,  they  should  be  large  enough  or  small  enough  that  they  can  be  measured  with 
reasonable  accuracy.  Two,  the  experimental  setup  should  be  fairly  straightforward  and 
simple.  Three,  the  analysis  should  not  be  overly  complex. 
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After  weighing  these  factors,  it  was  decided  that  relating  the  input  impedance  of  a 
dipoie  near  the  ground  to  its  permittivity  would  be  a  practical  experimental  technique.  The 
input  impedance  of  a  half-wave  dipole  can  be  expected  to  be  around  75  ohms,  well  within 
the  measurement  range  of  modem  RF  impedance  meters.  The  experiment  requires  only  a 
one-port  measurement  of  a  simple  antenna.  The  analysis  involved  is  manageable. 

The  method  chosen  for  theoretical  analysis  is  an  image  approach  developed  by 
Lindell  and  Alanen  [13,14,15,16].  They  replaced  the  ground  with  an  equivalent  image 
current  located  on  a  line  in  complex  space. 

In  Chapter  2,  Lindell's  image  current  is  rederived  for  the  special  case  of  a  horizontal 
electric  dipole  (HED)  over  a  ground  with  complex  permittivity  e  =  £r€0  and  magnetic 
permeability  p.  =  po-  It  is  then  applied  to  obtain  an  expression  for  the  input  impedance  of 
an  electric  dipole  of  arbitrary  length  over  ground. 

Numerical  results  are  given  in  Chapter  3,  convergence  of  integrals  is  shown,  and 
the  limitations  of  some  assumptions  will  be  addressed 

In  Chapter  4,  the  relationship  between  the  permittivity  of  the  ground  and  dipole 
impedance  is  discussed.  Graphs  are  shown  depicting  the  effects  of  dipole  height  and 
length  on  the  impedances. 

Finally,  Chapter  5  addresses  experimental  considerations. 
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2.  DERIVATION  OF  IMPEDANCE  EXPRESSIONS 


The  derivation  of  the  image  current  begins  with  the  two-dimensional  transform  in 
the  x-y  plane  of  Maxwell's  equations.  The  source  will  be  an  infinitesimal  HED  with 
current  magnitude  I,  length  L,  direction  x,  and  location  at  a  height  z  =  h  above  the  ground. 
Thus,  after  the  transform  of  Maxwell's  equations  with  respect  to  x  and  y,  we  have 


2xE'-jKxE  +jcopoH  =0 

z  >  0 

(2.1) 

2xH'  -  jKxH  -  jcoeoE  =  IL*5(z-h) 

z  >  0 

(2.2) 

2xE{  -jKxEj  +jcojioHi  =0 

z  >  0 

(2.3) 

2xHi  -jKxHi  -jcoeiEi  =0 

z  >  0 . 

(2.4) 

where  K  is  the  two-dimensional  transform  variable,  the  prime  indicates  the  derivative  with 
respect  to  z,  ei  »  er€o,  and  a  subscript  of  one  denotes  the  earth.  Note  that  the  magnetic 
permeability  is  assumed  to  be  the  same  as  that  of  free  space,  since  that  is  usually  the  case  in 
nature.  Equations  (2.1)  thru  (2.4)  can  be  uniquely  expressed  in  terms  of  the  transverse 
components  of  electric  and  magnetic  fields,  represented  here  by  e  and  h.  Therefore,  after 
some  manipulating  and  considering  only  the  transverse  fields,  we  are  left  with 


e'=jj-r(itk2-KK)xi-  h 

(2.5) 

h'as-jp(Itk2-KK)xie 

(2.6) 

=jj^*(Itki  -  KK)  x  4  •  h| 

(2.7) 

hi  =  -jj-*(  Itkj  -  KK)  x  2  ■  ei  . 


(2.8) 
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The  two-dimensional  idem  dyadic  is  It.  Also,  q  =  V  M-o/Eo.  q  i  =  Vdo/Eh  k  =  coV |i0£o. 
and  ki  =  cayj  (iod  .  For  the  sake  of  clarification,  dyadics  will  be  represented  by  double 
underlined  capital  letters  except  at  times  when  their  meaning  is  better  conveyed  by  writing 
two  vectors  next  to  each  other,  e.g.,  VV. 

Note  that  the  above  equations  are  of  the  same  form  as  the  telegrapher's  equations. 


U'(z)=-zdl(z)  +  u  (z) 
I'(z)  =  -  ydU(z)  +  i(z) 


(2.10) 


This  similarity  opens  the  possibility  of  representing  each  medium  by  an  equivalent 
distributed  impedance.  If  the  voltage  represents  e  and  the  current  h,  the  corresponding 
distributed  impedance  and  admittance  are  dyadics  given  by 


Zd  =  -jTl(Itk2-KK)xi  z  >  0 

Y^  =  ( Ijk2  -  KK)  x  t  z  >  0 

A  subscript  of  one  is  added  when  z<0.  The  product  of  Zd  and  Yd  is 

Zd- Yd  =  -  (k2  -  K2)  =  -32  It .  z>0  . 

Here,  3  is  the  z-component  of  the  wave  number,  k,  given  by 

3  =  Vk2  -  K2  ,  3l  =  Verk2  -  K2  =  V  (£r-l)k2  -  32 


(2.11) 

(2.12) 


(2.13) 


(2.14) 


and  K  is  the  magnitude  of  the  vector  K.  The  branches  of  3  and  3i  are  taken  such  that  their 
arguments  are  between  -n  and  0.  The  reason  for  this  choice  will  be  made  clear,  as  it  will 
insure  that  the  electric  field  approaches  zero  as  z— **>  in  Equation  (2.24). 

Now  that  the  equivalent  impedance  is  known,  the  next  step  is  to  find  a  suitable 
reflection  dyadic. 


R  =  (Z|  -Z)-(Zi  +  Z)-1 


(2.15) 
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The  inverse  of  a  two-dimensional  dyadic  is  the  same  as  that  of  its  associated  matrix. 


where  adj(A)  and  det(A)  are  the  adjoint  and  determinant  of  the  associated  matrix.  Applying 
(2.16)  to  the  inverse  dyadic  in  (2.15)  and  simplifying  the  result  yield  the  following 
reflection  dyadic: 


r  =  rTM  KK/K2  +  RTC  adj(KK)/K2 . 


(2.17) 


R™  and  RTE  are  the  familiar  reflection  coefficients  for  the  transverse  magnetic  and 
transverse  electric  cases.  They  are  given  by 

(2.18) 

P+Pi  p  -1  V  (ef-l)k2+{J2 

rtu  £rMl  £rp  -  V  (£r-l)k2+gi  (2  19) 

£rP+Pi  erP  +  V  (er-l)k2+p2 


The  general  form  of  this  reflection  dyadic  was  derived  by  Lindell  in  [13].  Here  RTE  has 
been  made  more  specific  by  setting  H=Ho-  The  adjoint  operator  and  the  denominator  K2 
make  the  reflection  dyadic  in  (2.17)  somewhat  awkward  to  work  with.  Applying  the 
dyadic  property 

adj(KK)  =  K2|jt  -  KK  ,  (2.20) 


R  can  be  written 


R  =  RTE  I,  +  Ro  KK/k2  , 

with  Ro  =  k2(R™  -  RTE)/K2  =  ?(^fr.I)k2  RTM  + 

(P+Pl)(erP+Pi)  Er  (  er+l 


(2.21) 

(2.22) 
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This  form  isolates  the  transform  variable  K  to  a  dyadic  in  only  one  of  the  terms.  Thus,  the 
inverse  transform  of  the  dyadic  pan  of  R  is  accomplished  by  simply  replacing  KK  with 

-VtVt,  where  Vt  is  the  x-y  component  of  the  vector  operator  V. 

Substituting  h'  from  (2.6)  into  (2.5),  a  Helmholtz-like  equation  is  obtained. 

e"  +  p2e  =  jco|i<>IL(  Itk2  -  KK)/k2  5(z-h) .  (2.23) 

From  this,  the  electric  field  is  found  to  be 

e  =  ^2- IL(  ^  -  k~2KK)  •  x  ■  ( ^  e"iP|2+h|  +  R  e~^z+h))  z  >  0  .  (2.24) 

Although  h  can  be  derived  similarly,  it  will  not  be  necessary  to  discuss  both  e  and  h.  Note 
that  the  first  term  above  is  that  portion  of  the  field  associated  with  a  dipole  in  free  space  and 
is  travelling  away  from  the  z  =  h  plane.  The  second  term  is  travelling  upward  and  can  be 
thought  of  as  a  wave  reflected  from  the  earth. 

This  leads  us  to  represent  the  electric  field  as  the  sum  of  the  well-known  free-space 
field,  Ei,  and  the  electric  field  scattered  from  the  ground,  E2. 

E  =  Ei+E2  (2.25) 

Ei  =  -jo)Ho  J  G  (lr+2l)-*I(x)  dx  ,  (2.26) 

where  G(r)  is  the  free-space  Green's  dyadic  function  ( I-VV)G(d),  and  G(d)=(e~jkd)/47td. 
Notice  that  taking  the  inverse  Fourier  transform  of  the  first  term  of  (2.24)  yields  only  the 
transverse  components  of  Ei.  Yet,  the  sole  difference  between  Ei  and  the  transverse  fields 
is  that  (2.24)  produces  the  dyadic  operator  ( It  +  VtVt)  instead  of  (  I-VV).  This  simple 

relationship  will  be  called  on  later  to  recover  the  z-component  of  E2  rather  than  carry 
longitudinal  components  through  the  derivation. 
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At  this  point,  instead  of  taking  the  inverse  Fourier  transform  of  (2.25)  (as  is 
normally  done),  the  reflection  dyadic  will  be  represented  as  the  Laplace  transform  of  a 
function: 

oo 

R(p)  =  J  F  (t)  e-P‘  dt .  (2.27) 

o 

After  inserting  (2.27)  in  (2.25),  changing  the  order  of  transformation,  and  taking  the 
inverse  Fourier  transform,  the  solution  of  E2  becomes 

00 

E2  =  -  jotyioIL  J  £(r+i(h-jt))-F(t)-x  dt  ,  (2.28) 

o 

which  implies  that  the  ground  can  be  replaced  by  an  equivalent  image  current 

OO 

E2(r)  =  -jcono  |  (^(r-z’i)-I(z')  dz’ 

I(z’)  =  -jIL  F(-j(z’+h»  *  (2.29) 

z'  =  jt  -  h  . 

What  remains  to  the  solution  is  to  find  an  expression  for  the  dyadic,  F(t).  If  we  solve 
Equation  (2.27)  for  F(t),  we  find  that  the  inverse  Laplace  transform  can  be  broken  into 
two  parts: 

F(t)  *  L-MR^JIt  -  L-1{R0}  KK/k2  .  (2.30) 

The  inverse  transform  of  can  be  found  in  the  Laplace  transform  tables  [18]  or  [19]  to 
be 

-2/;(taVeH) 

t 


L-l(RTE)  = 


(2.31) 
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-2y2(-jk(z'+h)Ve^I) 


z  +h 


(2.32) 


Thus,  the  contribution  of  the  first  term  of  R  in  (2.22)  to  the  image  current  is 

T  ,  ,  -2IL-/2(-jk(z'+h)Ve^I)  . 

W)  -  zVh  z'+h  ■  x. 


(2.33) 


On  inspection,  it  is  noted  that  if  z'  in  (2.32)  is  real  the  Bessel  function  will  have  a  complex 
argument  and,  therefore,  become  singular  as  z'  approaches  infinity.  In  fact,  to  insure 
convergence,  the  argument  of  the  Bessel  function  must  be  real.  Hence,  it  is  set  equal  to  the 
positive  real  variable  p: 

p  =  -  jk(z'+h)  .  (2.34) 


This  forces  z'  to  be  a  complex  number  given  by 


(2.35) 


Using  this  definition  of  z',  the  electric  field  due  to  this  image  current  can  be  written  as 


E2a  =  2jcopo  J  G  (d)  x  '^dp  .  (2.36) 

o 


The  branch  of  the  distance  function,  d,  is  defined  as  -n  <  arg(d)  <.  0  to  insure  that  the 
Green's  function  asymptotically  approaches  zero  allowing  the  integral  of  E2a  to  converge. 
With  this  definition  of  the  branch  of  d,  the  integral  in  (2.36)  is  fairly  well  behaved.  In  fact, 
it  decays  exponentially  after  p  is  large  enough  to  dominate  the  exponent  of  the  Green's 
function.  Therefore,  the  integration  can  be  performed  numerically  using  a  simple  computer 
routine.  The  numeric  evaluation  of  this  integral  will  be  discussed  in  further  detail  in  the 
next  chapter. 
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The  inverse  transform  of  the  second  term  of  the  reflection  dyadic,  Rq,  is  somewhat 
more  difficult  than  the  first  one.  Because  Ro  is  essentially  R™  plus  a  constant,  we  will 
begin  by  taking  the  inverse  transform  of  R™.  As  a  matter  of  convenience,  the  transform 
variable  (3  will  be  replaced  by  qk^Er-l  to  give 


r,T»  M 

6rq  -  Vq2+1 


(2.37) 


where  q  can  be  thought  of  as  a  dimensionless,  normalized  propagation  constant.  After 
rearranging  (2.37)  and  applying  the  Laplace  transform  property  of  differentiation,  R™  can 
be  expressed  in  the  form 


RTM=&d^+^(qW^T)2 

6r+l  q2-^2  q2-^2 


(2.33) 


Here,  y  is  defined  as  1  e^-1 .  The  first  term  above  corresponds  to  a  delta  function  in  the 


spatial  domain.  The  second  is  the  transform  of  the  hyperbolic  sine.  The  inverse  transform 
of  the  third  term  can  be  expressed  as  the  convolution  of  a  hyperbolic  sine  and  a  Bessel 
function. 


f(p)  =  L-MRTM} 

P 

=  5(p) - CrY sinh(Yp)  +  2ery  f sinh(y(p-s)) dp.  (2.39) 

6r+l  6r+l  J  p 


The  above  transforms  and  transform  properties  can  be  found  in  [17]  and  [18]. 

At  first  glance,  the  hyperbolic-sine  functions  make  f(p)  appear  to  diverge  as  p 
approaches  infinity.  However,  the  second  and  third  terms  tend  to  cancel  for  large  p, 
causing  f(p)  to  approach  zero  in  the  limit  as  p  approaches  infinity.  Yet,  these  terms  also 
create  significant  round-off  error  making  the  form  of  f(p)  in  (2.39)  impractical  for 
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numerical  computation.  Hence,  another  representation  of  f(p)  will  be  necessary.  Taking 


advantage  of  the  property 

oo 

J  f(p)e-YP  dp  =  E-^~ , 

0  £f+l 

(2.40) 

f(p)  can  be  written  as 

f(p)  =  $(p)  +  fE(p) 

£r+l 

(2.41) 

OO 

f£(p)=Je-rs(^.M)ds. 

(2.42) 

Here  we  introduce  the  function  fe(p)  as  the  nonsingular  part  of  f(p).  Note  also  that 
(£r+l)fE(p)/er  the  inverse  transform  of  Ro.  For  ease  of  numerical  calculation,  Lindell 
[16]  expressed  the  integral  (2.42)  as  a  couple  of  series  expansions. 


(2.43) 


fe(P)  —  2e^y3  e-1p  +  — ^~=sin(p+7t/4)  +  cos(p+Jt/4) 

EtpWyK 


15  /  23  9  4 \ . 

- r=-|“  TT  +  T~2  — 4  sm(p+it/4) 

*3pWST  32  2E,2  e? 


105  /  91  59  7  4  \  ... 

■““T - f=-  “  177  +  - 7 - 7  +  “T  COS(p+7C/4 

2£r3p9/2V^r  128  16E,2  er4  ef 


(2.44) 


The  series  (2.43)  is  an  exact  representation  of  fE(p).  However,  it  does  not  begin  to 
converge  until  2n  >  p.  Thus,  (2.44),  an  asymptotic  expansion  of  fE(p),  serves  as  a  better 
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means  of  calculating  it  for  large  p.  Reference  [15]  describes  a  method  that  can  be  used  to 
obtain  the  full  series  representation  of  (2.44).  It  was  found  by  trial  and  error  that  p  =  5  is  a 
fairly  good  point  to  make  the  transition  between  (2.43)  and  (2.44). 

Note  that  for  the  special  case  of  er  =  1,  fE(p)  becomes  fi(p)  =  -2J(p)/p.  This  is  the 
same  as  the  inverse  transform  of  R1^.  It  is  also  convenient  to  mention  at  this  time  that  fe(p) 
satisfies  the  differential  equation 


(2.45) 


This  is  seen  by  applying  the  transform  property  of  differentiation  in  conjunction  with  the 
inverse  transform  of  (2.37). 

The  dyadic  F(p)  can  now  be  expressed  in  terms  of  the  two  functions  fi(p)  and 

fe(p)- 

£(p)  =  fi(p)  I  -  —  fe(p)  VVk-2  .  (2.46) 

Er 


It  follows  that  the  scattered  electric  field  is 


E2  *  -jtopoIL  J  G^r-|r'-  ■  -ji —  z  jj  •  F^p)  •  *  dx  . 


Therefore,  the  image  current  is 


Ii  =  F(p).*IL 

=  ^fl(p)I  +  £^fc(p)VVk-2j  •  x  IL  . 


(2.47) 


(2.48) 


If  (2.48)  is  generalized,  to  allow  the  current  to  be  a  function  of  x,  the  image  becomes 


Ii  =  fi(p)I(x)x  - 


j(er»D 
£rk  yj  £| —  l 


fe(p)IXx)i  . 


(2.49) 


where  -jk'y  ej-1  dz  =  dp.  Thus,  the  image  current  associated  with  an  HED,  depicted  in 
Fig.  2. 1  (all  figures  appear  at  the  end  of  each  chapter),  consists  of  a  horizontal  and  a 
vertical  component  both  along  the  same  line  in  the  complex  z-plane.  As  would  be 
expected,  both  components  decay  to  zero  as  p  becomes  large. 

Finding  an  expression  for  the  input  impedance  of  the  HED  defined  in  Fig.  2.2 
using  the  image  current  above  is  a  straightforward  procedure  using  image  currents.  The 
input  impedance  is  defined  by  the  stationary  functional 


Z  = 


•  x  I(x)  dx  . 


(2.50) 


The  electric  field,  given  in  (2.26)  and  (2.47),  can  be  plugged  into  (2.50)  to  give 


Z  =  J~-jj  I(x)  *-G(lr-r'l)-xI(x')  dx’  dx 


jkrn 


J  I(x)4  G(lr-(r,-jp/Ver-i)0'L(P)  ^I(x  )  dp  dx' dx  .  (2.51) 


Equation  (2.51)  implies  that  the  impedance,  like  the  electric  field,  can  be  expressed  as  the 
sum  of  the  free-spacc  impedance  and  the  change  in  impedance  due  to  the  ground.  That  can 
be  done,  of  course,  only  if  the  current  distribution  is  not  significantly  changed  by  the 
ground.  The  validity  of  this  assumption  will  be  discussed  in  Chapter  3.  After  evaluating 
the  dyadics  and  applying  the  differential  equation  (2.45),  AZ  can  be  written  as 


AZ=M 

'oj 


f 


I(x)I(x')Ki(lx-x'l)-I  '(x)I'(x')f ■^Ke(lx-x'l)-^*--< 

Ur  £r+l 


)jdx' 


(D0)  px'dx,  (2.52) 


where 
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Kc(lx-x’l)  =  j  fe(p)G(D)  dp 
o 


(2.53) 


Ki(tx-x’l)  =  Jfi(p)G(D)  dp 
o 


(2.54) 


D(p)  =  A  /(x-x')2  +  f2h  - 


(2.55) 


k  V  £r~l 


D0  =  D(p=0)  =  V(x-x')2  +  4h2 


(2.56) 


If  a  sinusoidal  current  is  assumed  for  the  case  of  a  half-wave-dipole,  the  impedance  change 


AZ=jkn  J  cos(kx)cos(kx')Ki(lx-x’l) 
J-X/4 

-*/4 


-  sin(kx)sin(kx’)  p4Ce(lx-x'l)  -  ^-^-G(D0)1  dx'  dx  .  (2.57) 

Ur  er+l  J 


A  change  of  variables  of  0  =  klx-x'l  and  <}>  =  klx+x'l  allows  one  of  the  integrals  in  (2.57)  to 
be  done  analytically  leaving 


AZ  =  jkr\  J  (sin0  +  (jc-0)cos0)Ki(Ix-x'I) 


+  (sin0  +  (Jt— 0)cos0)  f—  Ke(tx-x'l)  -  E— G(D0)1  dx'  dx  .  (2.58) 

ler  er+l 


The  integrand  of  (2.58)  is  very  well  behaved,  and  is  not  difficult  to  integrate  numerically. 
For  a  dipole  of  arbitrary  length,  2H,  and  current  distribution,  sin(k(H~lzl)),  the 


impedance  change  is 
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kH 


' 

KIq 


n 


£r~  1 


AKi(0/k)-B  ■— Ke(lx-x’l)  -  — -G(D0) 
l  {&  £r+l 


d0 


jm. 

1 dl 


2kH 

r  (  r 

CKt(0/k)-D 


kH 


V 


-Ke(ix-x'l)-e-^-G(D0)) 

v£r  £r+l  ‘  )j 


d9 


(2.59) 


with 

Io  =  sin(kH)  (2.60) 

A  =  -(cos2kH+0sin2kH+l)sin0+(0+(2kH-20)cos2kH)cos0+sin(2kH-0) 

B  =  (cos2kH+l)sin0+(0+(2kH-20)cos2kH+sin2kH)cosO+sin(2kH-0) 

+  sin2kHcos(2kH-0) 

C  =  (2kH-20)  cos0  -  cos2kH  sin(2kH-0) 

D  =  (2kH-20)  cos0  +  cos2kH  sin(2kH-0) . 

Equation  (2.60)  shows  that  the  expressions  representing  the  coefficients  A-D  are  long. 
However,  they  are  much  more  efficient  than  the  alternative  which  is  the  double  integration 
of  (2.58). 

As  the  dipole  height  increases,  the  integral  (2.58)  can  be  approximated  by  an 
algebraic  expression.  For  large  kh,  D(p)  becomes  2h  -  jp/V^r-1  •  Hence,  by  applying  the 
integral  identity 


oo 


Jfi(p)  e-P/V^dp  = 


i-Vfr 

iW^ 


(2.61) 


The  new  Green's  function,  Ki,  can  be  simplified  to 


Kis- 


i-Vgf 


i+V* 


.-J2kh 


2nh 


kh  — »  oo 


(2.62) 
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Since  the  Green's  function  Ke  is  no  longer  a  function  of  0,  the  double  integration  for  the 
currents  can  be  performed  analytically.  The  resulting  expression  is  equivalent  to  the  image 
reflection  coefficient  approximation. 

e-j2kh  kh  ->  oo  (2.63) 


AZ~ 


-ih 


1-Ver 


2jtkh 


1+ 


Image  current  due  to  HED  over  an  infinite  half  space  of 
complex  permittivity,  e.  The  current  lies  in  the  complex 
z-plane  and  has  two  pans:  Iix  parallel  to  the  dipole,  and 
oriented  vertically. 
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3.  NUMERICAL  CONSIDERATIONS 


In  this  Chapter,  the  graphical  characteristics  of  some  of  the  equations  given  in 
Chapter  1  are  presented.  To  begin,  the  transcendental  function  fe(p)  is  studied.  This  is 
followed  by  a  discussion  of  the  integrand  of  Ke(8)  and  its  integration  to  infinity. 

The  function  fe(p),  represented  by  Equations  (2.43)  and  (2.44),  was  coded  as  a 
Fortran  subroutine  to  be  used  later  in  calculating  dipole  impedances.  As  a  test  of  the 
subroutine,  Figures  3.1  through  3.5  were  generated  and  compared  with  those  published  by 
Lindell  [15,16].  Figures  3.1  through  3.4  agree  quite  well  with  LindeU's.  These  graphs 
demonstrate  the  general  appearance  of  fe(p).  Notice  that  for  small  er,  the  functions  are 
much  larger  near  the  origin.  On  the  other  hand,  as  the  permittivity  increases,  the  magnitude 
of  fc(p)  is  more  evenly  distributed  along  the  p-axis.  In  fact,  when  er  =  50  (Figure  3.2), 
fe(p)  is  practically  a  line  with  a  small  downward  slope. 

Upon  close  inspection  of  Figures  3.1  through  3.5,  one  can  see  a  very  slight 
discontinuity  in  the  curves  at  p  =  5.  This  is  seen  best  in  Figure  3.2  when  Er  =  10.  It  is 
caused  by  the  inexactness  of  the  asymptotic  approximation  at  the  transition  between 
Equations  (2.43)  and  (2.44).  In  the  case  of  impedance  calculations,  this  bump  causes 
negligible  error.  However,  there  may  be  some  applications  that  require  higher  accuracy.  If 
this  is  the  case,  the  transition  between  the  two  representations  can  be  moved  to  a  higher  p, 
or  the  method  described  in  [15]  could  be  used  to  generate  the  next  few  terms  of  the 
asymptotic  expansion. 

Figures  3.5  through  3.7  show  how  the  magnitude  of  fe(p)  in  affected  by  the 
argument  of  er.  A  quick  comparison  of  Figure  3.5  and  LindeU's  [16]  results  (Figure  3.8) 
shows  a  discrepancy  for  large  arg(er).  To  determine  which  is  correct,  we  need  only  inspect 
Equation  (2.44),  the  first  few  terms  of  which  are 
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fe(p)  ~  -  2e^73  e_>p  + 


2V2 


erp 


3/2^ 


sin(p+jt/4)  + 


3(8-3cJ) 

2^p5«V2s 


XOS(p+7l/4) 


15  /  23  _9__ 

■&.W5T  "  31  a?"  < 


sin(p+rt/4) . 


(3.1) 


We  see  that  the  second  term  and  higher-order  terms  tend  to  zero  for  all  values  of  er  as  p 
becomes  large.  The  first  term  decays  with  an  attenuation  constant  x  =  y  =  er2-l  j  1 ,  t 


>  0.  If  we  take  the  extreme  case  of  the  permittivity  being  completely  imaginary,  we  see  that 
x  =  0  and  fe(p)  approaches  the  expression 


f£(p)~2er2Y3e-tt\ 


Re(er)  =  0  .  (3.2) 


In  fact,  for  this  case,  fe(p)  approaches  a  helical  form  as  p  -» <*».  Therefore,  Figure  3.5  is 
the  more  accurate  graph.  At  first  glance,  it  may  be  disturbing  that  the  image  currents  do  not 
decay  to  zero  as  However,  fe(p)  can  be  thought  of  as  a  distribution.  In  a 

distribution  sense,  the  limit  as  p  approaches  infinity  can  be  treated  as  zero.  This  follows 
from  the  Riemann-Lebesgue  lemma  [19]. 

The  program  used  to  generate  fe(p)  was  then  tested  in  the  two  integral  properties 

OO 

i  fe(P)  dp  =  — —  (3.3) 

o  e,+l 

c-p/V£TTdp  =  (3.4) 

lWer 

The  integrations  were  performed  with  an  adaptive  trapezoidal  routine  that  was  truncated 
when  the  integrand  became  negligible.  The  results  varied  with  the  step  size  and  the  point  of 
truncation.  It  is  apparent  from  the  results  of  this  test,  shown  in  Table  3. 1  (p.  23),  that  the 
computation  of  fe(p)  as  well  as  the  integration  are  being  performed  properly. 


Jfl(p) 

o 


The  image  currents  are  obtained  from  f£(p).  The  horizontal  component  of  the  image 
current  is  equal  to  f i(p),  while  the  vertical  current  is  proportional  to  the  derivative  of  fe(p). 
Some  plots  of  the  vertical  image  currents  for  different  values  of  er  are  given  in  Figures  3.9 
through  3.12.  Recall  that  in  impedance  calculations  the  vertical  image  current  need  not  be 
calculated  since  f£(p)  is  differentiated  a  second  time,  and  a  substitution  can  be  made  for 
f£(p)  via  the  differential  Equation  (2.45)  eliminating  all  terms  containing  derivatives  of 

fe(p)- 

The  next  step  is  the  integration  to  compute  the  Green’s  function,  Ke(0), 

oo 

Ke(lx-x’l)  =  f  fe(p)G(D)  dp  (3.5) 

o 

D(p)  =  ^j(x-x')2  +  |2h  -  .  (3.6) 


At  first  glance,  it  appears  the  integration  to  infinity  could  be  the  source  of  some  difficulty. 
When  p  becomes  large  enough  to  dominate  D(p)  in  Equation  (3.6),  we  see  that  G(D)  can  be 
approximated  by 


G(D)  =  y£r~i 


g-pA/vT 

4tcp 


p  large.  (3.7) 


Hence,  when  the  approximation  (3.7)  applies,  G(D)  decays  exponentially.  The  function 
f£(p)  also  diminishes  as  p  increases.  It  is  apparent  that  the  integration  can  be  truncated  after 
(3.7)  is  valid  and  allowed  to  decay  a  few  time  constants. 

For  er  of  small  arguments,  D(p)  moves  through  a  minimum  slightly  before  G(D) 


begins  to  decay,  causing  a  spike  in  the  integrand.  This  is  illustrated  in  Figures  3.13 
through  3.15.  As  the  normalized  height,  hA,  of  the  dipole  decreases,  the  spike  tends  to 
increase.  However,  while  the  height  is  greater  than  zero  the  spike  will  remain  bounded. 


Figures  3.13  through  3.20  are  examples  of  the  integrand.  The  first  three  are 
identical  to  those  published  by  Lindell  [14].  The  next  five  plots  (Figures  3.16-20)  were 
added  to  show  the  effect  of  a  change  on  the  parameters  8,  hA,  arg(er),  and  lerl.  Figure 
3.16  shows  that  an  increase  in  9  moves  the  peak  further  from  the  origin.  The 
corresponding  integrands  of  Ki(0)  are  plotted  in  Figure  3.17.  Note  that  the  peaks  in  these 
two  figures  occur  at  the  same  places  along  the  p-axis.  As  shown  in  Figure  3.18,  the  peak 
in  the  integrand  decreases  as  hA  increases.  Figure  3.19  shows  that  the  peak  is  sharp  only 
for  small  arguments  of  er.  As  seen  in  Figure  3.20,  increasing  the  magnitude  of  er  also 
moves  the  peak  further  from  the  origin  while  smoothing  it  out  at  the  same  time. 

It  is  clear  by  inspection  of  Figures  3. 1 3  through  3.20  that  the  integration  producing 
Ke(0)  and  Ki(0)  can  be  truncated  after  a  certain  point,  p0,  with  negligible  error.  In  these 
graphs  we  see  that  pQ  varies  as  0,  hA,  and  er  change.  Therefore,  a  valid  means  of 
detecting  the  truncation  point  is  needed.  One  solution  to  this  problem  is  to  test  segments  of 
the  integrand  to  determine  whether  that  part  of  the  integrand  is  more  or  less  than,  for 
example,  0.01%  of  the  integration  to  that  point.  This  method  was  used  with  the  adaptive 
integration  routine  mentioned  earlier  and  found  to  be  very  effective. 

Another  means  of  predicting  a  good  truncation  point  was  derived  by  equating  the 
magnitudes  of  two  approximations  for  D(p),  lD(small  p)l  and  lD(large  p)l.  This  leads  to 

I V B2-i-k2h2l  =  Ikh  -jpi fyj er-ll  .  (3.8) 

As  a  result,  the  p  satisfying  (3.8)  should  be  in  the  region  slightly  before  G(D)  decays 
exponentially.  It  can  be  approximated  by 

Pi  =  IV^I  (8kh+0) .  (3.9) 

If  we  allow  G(D)  to  decay  a  few  attenuation  constants  beyond  pi  and  add  10  to  it  as  an 
extra  precaution,  we  should  be  able  to  truncate  the  integration  when 
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-  IVeTil 


(8kh  +  0  +  3)  +10. 


(3.10) 


Using  (3.10)  in  place  of  infinity  as  a  limit  of  integration  proved  to  give  good 
results.  In  fact.  Figures  3.13  through  3.20  were  all  generated  using  the  aforementioned 
adaptive  trapezoidal  routine  with  (3. 10)  as  an  endpoint.  Note  that  the  plots  are  all  relatively 
smooth,  even  in  the  region  of  the  spike  (Fig.  3.15),  and  in  each  case  the  integrand  had 
become  small  enough  before  truncation  that  it  could  no  longer  be  read  on  the  graph. 

Equation  (3.10)  indicates  that  as  lerl,  h,  or  0  become  large  the  integral  for  the 
Green's  function,  Ke(9),  takes  longer  to  compute.  Thus,  direct  computation  of 
impedances  using  Equation  (2.51)  becomes  impractical  when  pi  is  large.  In  the  case  of 
large  h,  Equation  (2.63)  can  be  used.  No  asymptotic  approximations  were  derived  for 
large  lerl  or  large  0.  It  is  the  opinion  of  the  author  that  IQ(0)  can  be  represented  by  an 
infinite  summation  much  as  fe(p)  was.  Similar  integrals  were  found  in  [20]  and  [21]. 
Unfortunately,  no  solution  for  integrals  of  form  identical  to  IQfO)  could  be  found. 
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Table  3.1.  Percent  error  of  adaptive  trapezoidal  integration  routine  when 
applied  to  the  integral  identities  in  Equations  (3.3)  and  (3.4). 


REAL 

IMAGINARY 

%  ERROR 

%  ERROR 

EPSILON 

EPSILON 

EQ.  3.3 

EQ.  3.4 

2 

0 

0.0656 

0.1168 

10 

0 

0. 1048 

0.0012 

20 

0 

0.0835 

0 . 0009 

50 

0 

0.0758 

0.0013 

2 

-0.  5 

0.2832 

0 . 0009 

10 

-2 . 5 

0.0296 

0.0012 

20 

-5 

0.0560 

0 . 0009 

50 

-12.5 

0.0732 

0.0013 

2 

-1 

0.2464 

0.0006 

10 

-5 

0. 1117 

0.0012 

20 

-10 

0 . 0606 

0.0009 

50 

-25 

0.0647 

0.0012 

2 

-2 

0.1578 

0.0009 

10 

-10 

0.0366 

0.0012 

20 

-20 

0.0648 

0.0011 

50 

-50 

0.0705 

0.0012 

2 

-4 

0.0547 

0.0018 

10 

-20 

0.0515 

0.0010 

20 

-40 

0.0578 

0.0013 

50 

-100 

0.0183 

0.0012 

2 

-10 

0.0308 

0.0014 

10 

-50 

0.0166 

0.0012 

20 

-100 

0.0089 

0.0013 

50 

-250 

0.0074 

0.0008 

to 


Figure  3.3.  Ife(p)l  for  different  e. 
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Figure  3.11.  Magnitude  of  vertical  image  current  for  different  values  of  er. 


Magnitude  of  the  integrand 
er  =  8,  and  kh  =  0.00001 
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in  Figure  3.14 
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Figure  3.16.  Magnitude  of  integrand  of  Kc(0)  for  different  0,  with 
h/X.  =  0.02,  £f  “  10-jO. 


Figure  3.18.  Magnitude  of  integrand  of  Ke(0)  for  different  hA,  with 
£r  =  2-jO  and  0  =  3n/8. 


SO  0 
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Figure  3.19.  Magnitude  of  integrand  of  Kc(0)  for  different  arg(er),  with 
lerl  =  2,  0  =  3ti/8,  and  h/X  =  0.02. 
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4.  IMPEDANCE  CHARACTERISTICS 


This  Chapter  discusses  dipole  impedance  characteristics.  Because  the  case  of  the 
half-wave  dipole  is  often  discussed  in  the  literature,  it  will  be  considered  first.  Equation 
(2.58)  is  used  to  compute  the  effect  of  the  ground.  This  is  added  to  the  thin-wire  free- 
space  input  impedance,  73.1  +  j42.5  ohms.  The  integration  in  Equation  (2.58)  is 
performed  using  the  same  adaptive  trapezoidal  routine  as  was  used  in  the  calculations  of  the 
preceding  chapter. 

Figures  4.1  through  4.4  were  generated  to  be  compared  with  results  already 
published  in  [22],  [23],  and  [16],  respectively.  The  resistances  compare  well  with  those 
published.The  reactance,  on  the  other  hand,  does  not  agree  with  that  documented  by 
Lindell  [16],  The  reason  for  this  is  not  known.  As  one  would  expect,  these  figures  show 
that  the  ground  has  more  of  an  effect  on  the  input  impedance  of  antennas  at  lower  heights. 
Note  also  that  there  are  nodes  at  heights  slightly  before  integer  multiples  of  a  quarter-wave¬ 
length.  Hence,  heights  of  A78  or  less  will  result  in  the  largest  change  in  impedance  and  thus 
will  yield  the  highest  resolution  when  measuring  permittivity. 

As  a  study  of  the  accuracy  of  the  approximate  impedance  given  in  Equation  (2.63), 
Figures  4.5  through  4.7  were  generated.  Each  figure  shows  the  real  and  imaginary  pans  of 
the  input  impedances  calculated  using  both  Equations  (2.60)  and  (2.63).  The  RCM 
(reflection  coefficient  method)  gives  accurate  impedances  at  larger  heights.  As  the  relative 
height,  hA,  of  the  dipole  decreases,  the  RCM  impedance  becomes  less  accurate,  until  at 
low  heights,  it  no  longer  even  resembles  the  correct  input  impedance.  Although  the  RCM 
approximation  may  be  useful  in  many  applications,  it  is  clear  that  it  is  not  appropriate  for 
use  in  determining  ground  permittivity  since  the  approximation  seems  to  be  valid  only 
when  the  impedance  is  least  affected  by  the  permittivity. 
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Calculations  were  made  using  Equation  (2.59).  The  results  were  plotted  as  the  ratio 
of  change  in  impedance  to  free-space  input  impedance  as  a  function  of  2kH  in  Figures  4.8 
through  4.15.  The  first  four  contrast  impedances  for  different  heights,  while  the  last  four 
show  the  effect  of  changing  er.  In  the  plots,  2kH  =  tc  corresponds  to  a  half-wave  dipole 
and  2kH  =  2jc  to  a  full-wave  dipole.  The  free-space  values  were  determined  using  the 
theoretical  expression  given  in  [24].  The  first  four  of  these  figures  reinforce  the  idea  that 
the  impedance  change  in  general  increases  as  the  height  decreases  for  antennas  of  any 
length.  Notice  that  as  the  dipole  approaches  a  wave  length,  AR/R  and  A X/X  tend  toward 
zero.  Also,  the  ratio  of  the  imaginary  parts  approaches  zero  when  2kH  is  small.  When  the 
relative  changes  in  input  impedances  are  zero  regardless  of  ground  permittivity,  it  is 
difficult  to  infer  a  permittivity.  Therefore,  it  would  be  best  to  use  dipoles  with  2H  <,  3X/4. 
In  the  case  of  short  dipoles,  the  effect  of  the  half  space  can  only  be  seen  in  the  real  part  of 
the  impedance.  This  allows  only  one  degree  of  freedom  to  determine  both  the  real  and 
imaginary  parts  of  the  permittivity.  It  would  be  best  to  use  antenna  lengths  that  show 
variations  in  both  the  real  and  imaginary  parts  of  the  impedance.  Therefore,  to  obtain  good 
resolution  between  the  input  impedance  and  complex  permittivity,  the  length  of  the  dipole 
should  be  in  the  region  >74  £  2H  <;  3X/4. 

The  next  five  figures,  4.16  through  4.20,  are  contour  maps  showing  the 
relationship  between  ground  permittivity  and  input  impedance  at  different  heights  above  the 
ground.  After  a  quick  inspection  of  these  figures,  one  can  see  that  the  contours  are  larger  at 
lower  heights.  They  are  all  well  behaved.  Once  the  input  impedance  is  known,  the 
complex  permittivity  can  be  easily  determined  using  a  contour  similar  to  those  shown  here. 
Only  at  the  lower  heights  does  the  contour  begin  to  fold  onto  itself.  If  a  measurement  falls 
into  this  region,  another  measurement  can  be  taken  at  some  other  height  or  frequency  to 
isolate  the  correct  permittivity. 


Input  resistance  (R)  of  half-wave  dipole  as  a  function  of 
relative  height(h/X)  above  various  grounds.  This  compares 
well  with  graphs  published  by  Miller  [22]. 
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Figure  4.2.  Input  resistance  (R)  of  half-wave  dipole  as  a  function  of 

relative  height(hA)  above  a  ground  with  £f  =  6  -  j0.03.  This 
compares  well  with  graphs  published  by  Lindell  [23], 
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Figure  4.3.  Input  resistance  (R)  of  half-wave  dipole  as  a  function  of 

relative  height(hA)  above  various  grounds.  This  compares 
well  with  graphs  published  by  Lindell  [16]. 


too 
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Comparison  of  input  impedances  for  half-wave  dipoles 
calculated  using  the  RCM  approximation,  Equation  (2.63), 

and  direct  computation.  Equation  (2.58),  with  e,  =  10  -  K) 
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Comparison  of  input  impedances  for  half-wave  dipoles 
calculated  using  the  RCM  approximation.  Equation  (2.63), 

and  direct  computation.  Equation  (2.58),  with  er  =  10  -  jlO. 
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Comparison  of  input  impedances  for  half-wave  dipole: 
calculated  using  the  RCM  approximation.  Equation  (2.63) 

and  direct  computation.  Equation  (2.58),  with  er  =  10-jU)0. 


2.5 


53 


Figure  4.8.  Change  in  input  resistance  (AR)  relative  to  free-space  input 
resistance  (R)  at  different  heights  above  the  ground  for 
different  length  antennas  (2kH)  with  er  =  5  -  j50. 


space  input 
ground  for 


Change  in  input  impedance  (AR)  relative  to  free-space  input 
impedance  (R)  at  different  heights  above  the  ground  for 
different  length  antennas  (2kH)  with  er  =  5  -  j5. 


input  reactance  (AX)  relative  to  free- space  input 
X)  at  different  heights  above  the  ground  for 
lgth  antennas  (2kH)  with  er  =  5  -  j5. 
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Figure  4.12.  Change  in  input  resistance  (AR)  relative  to  free-space  input 
resistance  (R)  above  various  grounds  for  different  length 
antennas  (2kH)  with  hA  =  0.02. 


Figure  413-  Change  in  input  reactance  (AX)  relative 
reactance  (X)  above  various  grounds  f< 

antennas  (2kH)  with  hA  =  0.02. 
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resistance  (AR)  relative  to  tree-space  input 
bove  various  grounds  for  different  length 
vi th  hA  =  0.02. 
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Figure  4. 15.  Change  in  input  reactance  (AX)  relative  to  free-space  input 
reactance  (X)  above  various  grounds  for  different  length 

antennas  (2kH)  with  h/X  =  0.02. 
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Figure  4.18.  Contour  map  relating  half-wave  dipole  input  impedance 
(R+jX)  to  complex  ground  permittivity  (er)  with  hA  =  0.094. 


100 
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5.  EXPERIMENTAL  CONSIDERATIONS 


In  measuring  the  input  impedance  to  determine  ground  permittivity,  there  are  a 
number  of  sources  of  error.  These  include:  1)  impedance  measurement;  2)  frequency 
measurement;  3)  dipole  height;  4)  dipole  length;  5)  free-space  impedance;  6)  roughness  of 
the  ground;  7)  inhomogeneities  in  the  ground;  8)  change  in  the  current  distribution  of  the 
dipole  due  to  the  ground;  and  9)  electromagnetic  scattering  from  environment  other  than  the 
ground. 

The  first  two  sources  of  error  are  determined  by  the  measuring  instrument.  If  the 
proper  equipment  is  chosen,  these  errors  can  be  relatively  small. 

The  height  of  the  dipole  is  critical,  particularly  at  the  lower  heights.  Figures  4. 1 
through  4.4  showed  that  the  derivative  of  AZ  with  respect  to  the  relative  height,  (hA),  of 
the  dipole  is  large.  Hence,  even  a  small  error  in  the  height  could  cause  a  large  error  in  the 
permittivity.  This  source  of  inaccuracy  can  be  lessened  by  not  placing  the  antennas  too 
close  to  the  ground. 

An  error  in  the  measurement  of  the  antenna  length  has  a  similar  effect  as  an  error  in 
the  frequency.  It  changes  the  value  of  2kH.  An  error  of  this  son  may  or  may  not  be 
critical  depending  on  2kH  and  the  permittivity. 

Accurate  determination  of  the  free-space  impedance  can  also  be  a  problem.  To 
circumvent  this,  it  could  be  measured,  or  a  contour  map  relating  the  permittivity  to  the 
difference  between  input  impedances  at  two  different  heights  could  be  used. 

The  extent  of  the  effect  of  ground  roughness  on  these  measurements  is  difficult  to 
quantify.  Of  course,  it  is  more  important  that  the  ground  be  smooth  directly  under  the 
antenna  than  farther  away. 
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The  effect  of  ground  roughness  on  these  measurements  is  twofold.  It  decreases  the 
certainty  of  the  height  measurement  and  changes  the  scattering  properties  of  the  half  space. 
For  best  results,  the  ground  should  be  as  smooth  as  possible  with  little  or  no  vegetation. 
Of  course,  it  is  more  important  that  the  ground  be  smooth  directly  beneath  the  dipole  than 
farther  away. 

Inhomogeneities  in  the  ground  could  be  the  source  of  considerable  error.  Care 
must  be  taken  to  insure  there  are  no  pipes  or  cables  buried  near  the  test  area.  Natural 
inhomogenieties,  such  as  stratification,  can  also  distort  measurements.  If  stratification  does 
not  occur  before  one  skin  depth, 

/ 

5=0) 

its  effect  on  the  impedance  will  be  small.  A  possible  means  of  compensating  for 
stratification  would  be  to  take  measurements  at  a  number  of  different  heights  and 
frequencies  and  use  a  parameter  optimization  routine  along  with  the  theory  presented  in 
[25]  to  determine  the  permittivity  and  depth  of  each  layer  of  earth. 

In  Chapter  1,  it  was  assumed  that  the  dipole  current  distribution  is  sinusoidal  and 
that  the  presence  of  the  earth  has  no  effect  on  that  distribution.  Of  course,  neither  of  these 
assumptions  is  true.  However,  they  do  provide  a  good  approximation.  According  to 
Lindell  [16],  the  fields  scattered  from  the  earth  induce  negligible  change  in  the  dipole 
current  if  h»5a/2,  where  6a  is  defined  to  satisfy  the  equation 

ka  =  2k6a  e_1/k5a  .  (5.2) 

here  a  is  the  radius  of  the  dipole.  Figures  5. 1  and  5.2  are  plots  of  hA  =  for  different 
frequencies  and  radii.  The  regions  above  these  curves  correspond  to  h>5a/2.  Finally,  the 
measuring  equipment  and  other  objects  should  be  placed  as  far  from  the  antenna  as 
possible. 
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Of  an  the  possible  sources  of  etror,  only  the  current  Attribution  is  unique  to  the 

theory  presented  here.  Most  other  experimental  procedures  will  be  prone  to  the  same  or 
parallel  problems. 


Relationship  between  5^2X  and  log(aA).  The  relationship 
h/X  »  8j/2X  must  hold  true  for  impedance  calculations  to 


Relationship  between  8a/2X  and  frequency  for  various 

antenna  radii.  The  relationship  h/X  »  must  hold  true 
for  impedance  calculations  to  be  valid. 


6.  CONCLUDING  REMARKS 


The  theory  between  the  electric  properties  of  the  ground  to  the  input  impedance  of  a 
horizontal  electric  dipole  has  been  presented.  The  necessary  numerical  computations  were 
discussed.  Then  patterns  of  the  input  impedance  were  illustrated  with  respect  to  the  height 
and  length  of  the  dipole.  Practical  considerations  were  discussed  Finally,  recommenda¬ 
tions  were  made  for  experimental  methods  of  determining  the  ground  permittivity  using  this 
theory. 

Attempts  were  made  to  experimentally  validate  the  techniques  mentioned  herein. 
However,  these  attempts  were  unsuccessful.  The  primary  reasons  for  this  are  thought  to 
be  parasitic  reactances  at  the  feed  point  and  scattering  from  equipment  used  in  the 
experiment. 

Despite  the  shortcomings  of  the  present  study,  it  is  the  opinion  of  the  author  that  the 
input  impedance  of  a  dipole  should  be  a  good — although  perhaps  not  the  best — means  of 
determining  ground  permittivity.  One  possible  improvement  to  this  technique  would  be  to 
use  a  slot  antenna  with  a  more  directional  radiation  pattern.  This  would  increase  the 
dependence  of  the  input  impedance  on  the  properties  of  the  earth  while  decreasing  its 
sensitivity  to  the  measuring  equipment  nearby.  The  parasitic  reactances  would  be 
decreased. 

An  improvement  that  could  be  made  to  the  theory  would  be  the  replacement  of  the 
integral  to  infinity  in  Equations  (2.53)  and  (2.54)  by  one  or  two  series  representations. 
This  would  decrease  the  computation  time  enough  to  make  it  practical  to  iteratively 
determine  the  current  distribution  on  the  dipole  before  using  the  stationary  functional 
(Equation  (2.50))  to  compute  the  input  impedances. 
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